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ABSTRACT: We analyse the quantum evolution of a particle moving in a potential 
in interaction with an environment of harmonic oscillators in a thermal state, using the 
quantum state diffusion (QSD) picture of Gisin and Percival. The QSD picture exploits a 
mathematical connection between the usual Markovian master equation for the evolution 
of the density operator and a class of stochastic non-linear Schrodinger equations (Ito 
equations) for a pure state \ifi), and appears to supply a good description of individual 
systems and processes. We find approximate stationary solutions to the Ito equation (exact, 
for the case of quadratic potentials). The solutions are Gaussians, localized around a point 
in phase space undergoing classical Brownian motion. We show, for quadratic potentials, 
that every initial state approaches these stationary solutions in the long time limit. We 
recover the density operator corresponding to these solutions, and thus show, for this 
particular model, that the QSD picture effectively supplies a prescription for approximately 
diagonalizing the density operator in a basis of phase space localized states. We show that 
the rate of localization is related to the decoherence time, and also to the timescale on which 
thermal and quantum fluctuations become comparable. We use these results to exemplify 
the general connection between the QSD picture and the decoherent histories approach to 
quantum mechanics, discussed previously by Diosi, Gisin, Halliwell and Percival. 



1. INTRODUCTION 



One of the basic premises of quantum theory is that the quantum state of a genuinely 
closed and isolated system evolves according to the Schodinger equation. Although some 
systems of interest are approximately closed and isolated, most of the systems we encounter 
are not, as a result of either purposeful intervention by measuring devices, or unavoidable 
interaction with the immediate environment. Such systems are said to be open, and are 
often studied in quantum optics [1], quantum measurement theory [2,3], and in connection 
with decoherence and emergent classicality [4,5,6]. 

An open quantum system is in essence a distinguished subsystem of a large, closed 
and isolated system in which there is a natural division into subsystem and environment. 
Although such divisions of the world cannot be explicitly identified in general, they do 
exist in a wide variety of situations of both experimental and theoretical interest. For 
example, in quantum optics, the distinguished subsystem is an atom or small collection of 
atoms, and the environment is the electromagnetic fields in interaction with it. We will 
in this paper be primarily concerned with that paradigm of open quantum systems, the 
quantum Brownian motion model, which consists of a large particle coupled to a bath of 
harmonic oscillators in a thermal state [7,8]. 

If the state of the total quantum system is described by a density operator Pt t a l then 
the state p of the subsystem is obtained by tracing it over the environment. An evolution 
equation for p (a master equation) may then be derived. This is in principle obtained 
quite simply by tracing the unitary evolution equation for ptotal over the environment. In 
practice, this is hard to carry out with any degree of generality, and has been carried out 
in detail only in specific examples (see Ref.[9] for example). As an alternative, one can ask 
for the most general evolution equation for p that preserves density operator properties: 
hermiticity, unit trace, and positivity. These conditions alone do not allow one to say very 
much about the form of the equation, but if one makes the additional assumption that the 
evolution is Markovian, then the master equation must take the Lindblad form [10], 



Here, H is the Hamiltonian of the open system in the absence of the environment (some- 
times modified by terms depending on the Lj) and the n operators Lj model the effects of 
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(1.1) 



the environment. For example, in the quantum Brownian motion model, there is a single 
non-hermitian L which is a linear combination of position and momentum operators. The 
Markovian assumption is not always valid, but is known to be a good approximation for 
a wide variety of interesting physical situations, e.g., for the quantum Brownian motion 
model in a high temperature environment. 

Density operators evolving according to a master equation (not always of precisely the 
above form) have been the subject of a number of studies concerned with decoherence and 
the emergence of classical behaviour [6,9,11,12,13,14,15,16]. In particular models, it has 
been shown that the density operator can become approximately diagonal in some basis 
(sometimes more than one), indicating that interference between the states in that basis is 
destroyed. This suggests that one has some right to regard the dynamical variables corre- 
sponding to the diagonalizing basis as "definite" . One may then ask for the probabilities 
of successive values of these variables, and whether those probabilities are peaked about 
approximately classical evolution. 

This approach to emergent classicality has considerable intuitive appeal, but there at 
least two ways in which it could be made more precise. 

First of all, the notion of diagonality of the density operator is too vague. In the 
quantum Brownian motion model, for example, one expects both position and momentum 
to become reasonably definite. The argument as to how this comes about often goes 
as follows [13]: The coupling to the environment is typically through position, and the 
density operator tends to become approximately diagonal in position very quickly. On 
longer timescales, the Hamiltonian part of the evolution begins to contribute, and the 
basis of diagonalization is rotated in phase space. As a result of this interplay between 
the Hamiltonian and the interaction with the environment, the density operator therefore 
becomes approximately diagonal in a basis of states that are localized in phase space, such 
as coherent states. It is, however, difficult to see this precisely and with any degree of 
generality. What is required is an explicit way of exhibiting the diagonality in phase space 
localized states. 

Second, the way in which one attempts to see the emergence of classical behaviour 
for the variables which have become definite is to consider the evolution of states initially 
localized in phase space. Such states will tend to follow approximately classical trajectories 
in phase space, with spreading due to quantum and environmentally-induced fluctuations. 
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The emergence of approximately classical trajectories is, however, much harder to see for 
arbitrary initial states, e.g., for superpositions of phase space localized states. The problem, 
in essence, is that the density operator does not in general correspond to the behaviour of 
an individual trajectory, but to an ensemble. The density operator for an arbitrary initial 
state will be very spread out in phase space, and it is not at all clear that it corresponds 
to the intuitive expectation of a statistical mixture of classical trajectories. Clearly what 
would be very useful in this context is an alternative description of open systems that 
could give a clearer physical picture of the behaviour of an individual system, rather than 
ensembles. 

A recently developed picture of open systems that may be the required precision tool 
is the quantum state diffusion picture, introduced by Gisin and Percival [17,18,19]. In this 
picture, the density operator p satisfying (1.1) is regarded as a mean over a distribution 
of pure state density operators, 

p = M\j,)Q>\ (1.2) 

where M denotes the mean (defined below) , with the pure states evolving according to the 
non-linear stochastic Langevin-Ito equation, 

j 

+ Y^{L j -(L j ))\i>)d£ j (t) (1.3) 
j 

for the normalized state vector Here, the d£j are independent complex differential 

random variables representing a complex Wiener process. Their linear and quadratic 
means are, 

M[d£ j dZt] = 6 jk dt, M[dZ j d£ k ] = 0, M[dtj] = (1.4) 

The master equation (1.1) is invariant under unitary transformations of the Lindblad 
operator, Lj^J^k ^jk^ki where Uj k are the components of a unitary matrix [17]. Physics 
therefore corresponds to the equivalence class of master equations equivalent under these 
transformations. Correspondingly, the Ito equation (1.3) is invariant under the same uni- 
tary transformations on the Lj's, supplemented by similar transformations on the noise 
terms, and thus there is an equivalence class of Ito equations also. 

The precise mathematical relation between (1.3) and (1.1) is that the class of Ito 
equations (1.3) is equivalent to the class of master equations (1.1). Indeed, this connection 
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supplies an alternative algorithm for numerical solution [19]. However, the strength of this 
picture is that solutions to (1.3) supply an intuitively appealing picture of the expected 
behavior of individual systems, and have been seen to correspond very closely to individual 
runs of actual laboratory experiments in quantum optics [20]. 

The connection between (1.1) and (1.3) is closely analagous to the connection be- 
tween the Fokker-Planck equation and the Langevin equation in the classical description 
of Brownian motion. There, one has two completely equivalent mathematical descriptions 
with very different pictures. The Fokker-Planck equation describes an ensemble of systems 
evolving deterministically, whilst the Langevin equation describes an individual system 
evolving stochastically. 

The quantum state diffusion picture has much in common mathematically with a vari- 
ety of recent attempts to modify quantum mechanics at a fundamental level [21,22,23,24]. 
In such attempts, equations of the form (1.3), or similar, are proposed. The difference be- 
tween QSD and such alternative formulations, is that QSD is regarded as a phenomenologi- 
cal picture, appropriate only under certain conditions, whereas the alternative formulation 
are taken to be fundamental. Eq.(l.l) and (1.3) also arise in descriptions of continuous 
measurement in standard quantum mechanics [25,26]. This paper is primarily concerned 
with the mathematical properties of Eq.(1.3), hence the results will be applicable to all of 
these situations. 

Solutions to the Ito equation often have the feature that they settle down to solutions 
of rather simple behaviour after a period of time. This general pattern of behaviour is 
indicated by numerical solutions [19], along with a number of localization theorems, which 
show that, unlike evolution according to the master equation, the dispersion of certain 
operators decreases as time evolves [18,27]. That is, certain types of variables become 
more "definite" as time evolves. 

A particularly useful example for our purposes was given by Diosi [28] , who considered 
the Ito equation (1.3) with L = ax and H = /2m. (This is the quantum Brownian 
motion model for the free particle neglecting dissipation). He showed that there exist 
stationary solutions \^ pq ) to the Ito equation that consist of Gaussian wave packets tightly 
concentrated about a point in phase space evolving according to the stochastic equations 
of classical Brownian motion. This is a particularly appealing result. The solutions to the 
Ito equation correspond very closely to macroscopic observations of an individual particle 
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interacting with an environment. 

Given a set of localized phase space solutions \^pq), such as the Diosi solution, the 
density operator may be reconstructed via (1.2). This, it may be shown, may be written 
explicitly as 

p = J dpdq f{p,q,t)\y pq ){y pq \ (1.5) 

where f(p, q,t) is a non-negative, normalized solution to the Fokker-Planck equation cor- 
responding to the Langevin equation describing the Brownian motion of the centre of the 
stationary solutions. 

The crucial point, now, is that the representation (1.5) of the density operator provides 
the desired improvements of the density operator programme described above. Firstly, the 
stationary states \^ pq ) are approximately orthogonal (for sufficiently distinct values of their 
centres, p,q). Eq.(1.5) therefore shows explicitly how the density operator may achieve a 
form in which it is approximately diagonal in a set of phase space localized states. Secondly, 
each diagonal element corresponds to an individual classical trajectory (with noise) . This 
means that the density operator might reasonably be interpreted as corresponding to a 
statistical mixture of classical trajectories. 

The object of this paper is demonstrate the above statements in detail, for systems 
more general than the case considered by Diosi. We will consider an open system consisting 
of a particle moving in a potential V(x), coupled to an environment described by Lindblad 
operators in (1.1) which are a linear combination of position and momentum operators. 
The detailed description of the model is given in Section II. 

We shall show, in Section III, that the Ito equation (1.3) has stationary solutions con- 
sisting of Gaussian wavepackets concentrated about points in phase space which undergo 
classical Brownian motion. These solutions are exact for quadratic V(x). The solutions 
for general potentials V(x) are approximate, and are valid as long as the localization width 
is much smaller than the length scale on which the potential varies. 

We shall then show, in Section IV, that every initial state tends towards one of the 
stationary solutions, for linear systems. In Section V, we consider the rate of localization, 
and show that it is related to the decoherence time, and also to the timescale on which 
thermal and quantum fluctuations become comparable. 

In Section VI, we construct the density operator of the form (1.5) explicitly, and discuss 
its properties. 
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Arguably the most comprehensive and fundamental approach to the problem of emer- 
gent classicality in quantum theory is the decoherent histories approach [4,29,30,31]. In 
fact, in Ref.[32], it was argued that there is a close connection between the quantum state 
diffusion picture and the decoherent histories approach. In Section VII, we use the above 
results to elaborate on this connection. 

We summarize and conclude in Section VIII. 



2. THE MODEL 



In this paper, we are concerned with systems described by a master equation of the 
form (1.1) with a single non-hermitian Lindblad operator linear in x and p 

L = ax + ibp (2.1) 

where a and b are real constants. The unitary transformations under which the master 
equation is invariant reduce to a simple phase invariance, L—>e™L. What follows therefore 
applies also to L's of the form (2.1) multiplied by a phase. This form of L is sufficient to 
describe the quantum Brownian motion model (see below), but also includes the cases in 
which L is taken to be a creation or annihilation operator. 

The operator H in (1.1) is taken to be 

H= ^— + V(x) + c{x,p} = H + c{x,p} (2.2) 

where c is a real constant. The master equation may then be written, 

i 1 1 o 1 o 

p = --[H + (c- -hab){x,p},p] -iab[x,{p,p}} - -a [x,[x,p]]--b \p,\p,p]} (2.3a) 

or alternatively, 

% 1 1 2 1 2 

p = --[Ho+(c + -hab){x,p},p]+iab\p,{p,x}]--a [x,[x,p]]--b \p,]p,p]] (2.36) 

Hereafter, we take c = \abTi. This ensures that the Ehrenfest type result, Tr(pp) = 
^Tr(£p), holds. 



The corresponding Ito equation is 

|#> = ~ (h + ±hab{x,p}J \1>)dt 

- - (a 2 (x - (x)) 2 + b 2 (p - (p)) 2 + 2iab(x - (x)p) - Tiab^j \ip)dt 

+ (a(x - (x)) + ib(p - (p))) |VR (2.4) 

We are particularly interested in the quantum Brownian motion model, for which the 
Lindblad operator is as above, but with 

a = (2£>)-3, b=(2D)*l, c=]p (2.5) 

Here, D = % / (8m^kT) , where 7 is the dissipation of the environment and T is its tem- 
perature. The master equation in this particular case may then be written, 

P = -j[Hq,p] - ^l[xApip}\ ~ 2il ^ fcT [x, [x,p\] - gi J fcr [p, \p,p\] (2.6) 

This does not, in fact, completely agree with the master equation given in a number of 
previous papers on quantum Brownian motion. In particular, the master equation given by 
Caldeira and Leggett [7] does not involve the term [p, [p, p]]. This difference is due to the 
fact the above master equation, by design, respects the positivity of the density operator, 
whilst the Caldeira-Leggett equation is known to violate it on short time scales [33] . This 
difference is not important, since we expect the Markovian approximation to hold only 
for high temperatures, and in this case the extra term is negligible since its coefficient is 
proportional to T _1 . (See Ref.[34] for further discussion, and also Ref.[9] for the derivation 
of exact master equations). 

Some information on the behaviour of the solutions to the Ito equation may be obtained 
by computing the time evolution of the moments of x and p, and this will be useful in the 
following sections. For any operator G, the time evolution of its expectation value in the 
state \ip) is given by 

d(G) = (ip\G\di{;) + {d*l>\G\*l}) + {d%l)\G\d%l)) 

= % -([H, G])dt - i(Lt[L, G] + [G, tf]L)dt 

+ a(G^,L)d£ + a(L,G)dC (2.7) 
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Here, following Percival [27], we have introduced the notation 

a(B, C) = ((flt _ {B )*)(C - <C»> = (Btc> - (B)*(C) (2.8) 

for the correlation between two operators £?, C in the state 

Setting G equal to p and x in this equation we obtain the Langevin equations 

d(x) = —dt + a(x, L)d£ + a(L, x)d£* (2.9) 
d(p) = -{V'(x))dt-2hab{p)dt + (j(p,L)dZ + (j(L,p)dC (2.10) 

With the choice of parameters (2.5), and for quadratic potentials, these equations describe 
classical Brownian motion. For more general potentials, this is true only if the state is 
sufficiently well- localized in x for the approximation (V'(x)) ~ V'((x)) to be valid (see 
below) . 

It is also of interest to compute the mean of time evolution of higher moments of x 
and ]3, and these may again be computed using (2.7). One finds, 

M d[^xf = 2R + 2hab{Ax) 2 + 2b 2 (^- R 2 ] - 2a 2 (Axf (2.11) 



dt m \ 4 

dt I 2 



= _ 2 (\($v'{x) + V'{x)p) - (p)(V'(x)) 



2hab(Ap) 2 + 2a 2 I - R 2 \ - 2b 2 {Apf (2.12) 



M^ = -((xV'(x))-(x)(V'(x))) + { -^ 

- 2a 2 R(Ax) 2 - 2b 2 R(Ap) 2 (2.13) 

Here, R is the symmetrized correlation between p and £, 

R = 1 {a{x, p) + a{p, x)) = a{p, x) + y = a{x, p) - y (2.14) 

Also, (Ax) 4 denotes ((x — (x)) 2 } 2 , and similary for (Ap) 4 . 

To handle general potentials is too difficult except in special cases, so approximations 
are required. Under Schrodinger evolution in ordinary quantum mechanics in a wide vari- 
ety of potentials, there exist approximate solutions consisting of localized Gaussian wave 
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packets concentrated about a classical path [35]. These solutions are possible because a 
sufficiently localized packet will only "notice" the quadratic approximation to the potential 
in the neighbourhood of the wavepacket's centre. The solution breaks down after a period 
of time, however, as a result of spreading of the wavepacket. 

Similar types of solution to the Ito equation (2.4) are possible, as we shall see in the 
next section. These have the advantage that wavepackets tend to localize with time, rather 
than spread. We may therefore justifiably approximate the potential-dependent terms in 
(2.12) and (2.13) by their expansions about the mean values of x and p. 

To see this more explicitly, and to assist the estimation of the validity of the approxi- 
mation, introduce the notation, x = (x), p = (p), and then write the potential as, 

V(x) = V(x) + (x- x)V'{x) + -(x- x) 2 V"{x) + W(x, x) (2.15) 

where 

W(x, x) = hx- xfv"\x) + ^(x- xfV^ + ■■■ (2.16) 
6 24 

Then the potential-dependent terms in (2.10), (2.12) and (2.13) become, 

(V'(x)) = V'{x) + (W'(x)) (2.17) 
(xV'(x)} - (x)(V'(x)) = (Ax) 2 V"(x) + ((x - x)W'(x)} (2.18) 

and 

\{pV'{x) + V'(x)p) - (p)(V'(x)) = RV"{x) + ^(pW'(x) + W'(x)p) - (j>)(W(x)) (2.19) 

The quadratic appproximation to the potential will therefore be valid when the terms 
involving W may be neglected in the above expressions. This will generally depend on the 
particular state. 

Taking the first few terms in the Taylor expansion of W, Eq.(2.17) for example, implies 
that the higher order terms may be neglected if 

V'{x) » l(Axf V"'{x) (2.20) 

This is clearly the condition that the width of the state is much less than the length scale 
on which the potential varies, as one would intuitively expect. The higher order terms in 
(2.18) and (2.19) also may be neglected if essentially the same type of condition holds. 
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3. STATIONARY SOLUTIONS TO 
THE LANGEVIN-ITO EQUATION 



We now show how to find stationary solutions to the Langevin-Ito equation, (2.4). It 
may be written 

\dij>) = u\ip)dt + v\ip)d£ (3.1) 

where 

i 1 

u = — + 7;foab + iab ((x)p — (p)x) 

-\a\x-(x)f-\b\p-(p)f (3.2) 
v = L-(L) (3.3) 

It is then convenient to rewrite the Ito equation in the exponential form 

|^) + \dtj}) = exp (udt + vd£) \ip) (3.4) 

The Diosi stationary solution has the feature that under time evolution, its shape is 
preserved and the only things that change are (ft) and (x) (and possibly a phase) [28]. Our 
approach to the search for stationary solutions to our more general equation is to require 
that the solution have this property. We therefore look for solutions to (3.1) satisfying the 
condition, 

m + |#> = exp (jrxd(p) - % -pd(x) + % -d(j^ (3.5) 

This is the statement that the state at time t + dt differs from the state at time t by nothing 
more than a phase, and a shift of (p) and (x) along the classical Brownian path described 
by (2.9), (2.10). Clearly (3.5) will be satisfied for any states of the form 

|V>) = exp {^(p) - ^P(x)^j \x) (3-6) 

where |x) is an arbitrary fiducial state. These are generalized coherent states [36]. 
We will solve (3.4) and (3.5) by first combining them to yield 



exp (udt + vd£) \ip) = exp ( —xd(p) — —pd(x) + —dtp ) \ip) (3-7) 

lb IV ft 
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and later confirm that the solution satisfies (3.5). 

Taking the operator on the right-hand side of (3.7) over to the left-hand side, and 
combining the exponentials using the Baker-Campbell-Hausdorff formula, one obtains, 

——xd{p) + —pd(x) — —dtp + udt + vd£ 

ft ft It 

~[x,v] d(p)dC + ^-[p,v] d(x)d£) |V) = IV) (3-8) 

Inserting the explicit expressions for d{p), d(x), u and v, and writing d(f> = fodt + <f>id£ + 
4>\d^,* (where 4>q is real), this equation becomes 

exp (Adt + Bdi + CdC) IV) = IV) (3-9) 

where 

A = u+ l - ({V'(x)) + 2hab(p)) x + ^p + 1<j(L, L) - U (3.10) 
a ami n 

B= l T (-<r(p, L)x + a(x, L)p - 0i) + L - (L) (3.11) 
n 

C= ^(-a(L,p)x + a(L,x)p-(/>l) (3.12) 

Expanding the exponential in (3.9), it follows that the state must obey the three equations, 

i|V) = (3.13) 
B|V)=0 (3.14) 
C\1>) = (3.15) 

Eqs.(3.14) and (3.15) will be satisfied if 

01 = a(x, L) (p) — a(p, L) (x) (3.16) 

and if the wave function is 

(a#) = Nexp (-/3(x - (x)) 2 + j-(p)x^j (3.17) 
for some constant (3, to be determined. The solution satisfies, 

(a#) + (x\dif>) = N exp (^-(3{x - (x) - d(x)) 2 + '-{(p) + d(p))x^j (3.18) 
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This is clearly a generalized coherent state, and thus satisfies Eq.(3.5). 

An equation for f3 may be obtained by inserting (3.17) in (3.13). One obtains the 
purely algebraic equation 



4 \b 2 + ^rj U 2 f3 2 + 4habf3 - (a 2 + ^V"({x))j = (3.19) 

where we have neglected terms higher than quadratic in the potential, as described in the 
previous section. 

It is of course possible to write down the explicit solution for /3, but it will generally 
be more useful in what follows to proceed differently. We have the uncertainty relation, 

h 2 

(Ax) 2 (Ap) 2 -R 2 >-j (3.20) 

with equality if and only if the state is of the form (3.17) [37]. Let us denote the values of 
the variances and correlation of the stationary state (3.17) by a 2 , a 2 and Rq. Then 



a 2 a 2 - R 2 = ^ (3.21) 



and 

fl-C-y (3.22) 

Since, from (3.19), (3 is a constant (to the extent that the approximation (2.20) holds) the 
stationary values of the variances and correlation must be those for which the right-hand 
sides of (2.11)-(2.13) vanish. That is, 

^° + haba 2 + b 2 (^-Rl]- a 2 a A x = (3.23) 



m J \ 1 

■ii,=.\t> *„*._2 , „2 ( h />:» \ ,,2 J 



-V"(x)Rq - habal + a 1 I ^— - R\ \ - b z a^ = (3.24) 



2 

-o%V"(x) + °^-- 2a 2 R a 2 - 2b 2 R§a 2 v = (3.25) 

These will be the most useful equations to work with in the following section. 

To see the complete solution in a particular case, let V(x) = and 6 = 0. The solution 
for (3 is then, 

i 



13 



where we have chosen the square root so that Re(3 > 0, for normalizability of the state. It 
follows that 

2 / Ti \ 2 9 ( h 2 ma 2 \ 2 Ti , 



This a close to minimal uncertainty state, since it satisfies, 

h 

V2 



&p&x = — (3.28) 



The solution (3.26)-(3.28) is very similar to the solution obtained by Diosi [28], but 
differs by some simple numerical factors, e.g., he obtained 

(°l)diosi = ' ( 3 - 29 ) 

This difference is due to the fact that Diosi used an Ito equation with a single real Wiener 
process, whereas the Wiener process used here is complex. 

The Diosi solution is also discussed in Ref.[38]. Some stationary solutions to (1.3) for 
the harmonic oscillator have also been found for by Salama and Gisin [39] , but their choice 
of Lindblad operators differs from that used here. 

Approximate stationary solutions to the Ito equation (2.4), for general potentials, are 
currently being studied by Brun et al. [40]. 



4. A LOCALIZATION THEOREM 



We now show that all solutions to the Ito equation tend towards the stationary solution 
in the long-time limit. The demonstration applies primarily to the case of linear systems, 
but we will work with a general potential in what follows, saving until the end the issue of 
the extent to which that case is properly covered here. 

We have shown that there is a two-parameter family of stationary solutions, parametrized 
by their centres (x), (p). To prove that all solutions tend to a stationary solution, we will 
exploit the fact that the stationary solutions are uniquely characterized by the statement 
that they are the eigenf unctions of the operator 

A = p-2ih(3x (4.1) 
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where (3 is the solution to Eq.(3.19). This means that the stationary solutions are uniquely 
defined by the statement that (AA) 2 = 0. We shall prove the desired result by showing 
that (AA) 2 tends to zero, in the mean. 

A number of "localization theorems" , showing that the dispersion of certain operators 
decreases with time, in the mean, have been proved by Gisin and Percival [18] and by Per- 
cival [27]. None of these results is applicable to the present case because their assumptions 
are too restrictive. They assume, for example, that the Hamiltonian is zero (or negligible), 
or that the Lindblad operators commute with the Hamiltonian. In brief, they assume that 
the Hamiltonian plays no significant role. An important feature of the case considered in 
this paper is that the stationary solutions are possible as a result of a balance between the 
wavepacket spreading induced by the Hamiltonian and the localizing effect of the Lindblad 
operators, and hence the role of the Hamiltonian cannot be ignored. An argument for the 
local stability of the stationary solution in the free particle case with 6 = was given by 
Diosi [28], but this proves nothing about arbitrary initial states. 

Returning to the case at hand, we have 



The rate of change of (AA) 2 in the mean, Md(AA) 2 , is then easily computed from 
Eqs.(2.11)-(2.13). It is convenient to write 



(AA) 2 = a (A, A) 



(Apf + 4h 2 \f3\ 2 (Ax) 2 - 2ih(f3 + P*)R - h 2 (/3 + (3*) 



(4.2) 



(Ax) 2 
(Ap) 2 



a 2 x (l + X) 
a 2 p (l + Y) 
R (l + Z) 



(4.3) 
(4.4) 
(4.5) 



R = 



hence the stationary solution is X = Y = Z = 0. One then obtains, 



M 



d(AA) 2 



= c x X + c 2 Y + c 3 Z 



dt 




+ Aa 2 Rl XZ + Ab 2 ^Rl YZ 



(4.6) 



x 
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where 



Cl = -h 2 a 2 + 2habal + 2R Q V"(x) (4.7) 
c 2 = -lUahal - 2 Jh°\_ h 2 b 2 ^ (4.8) 

c,= 2 -^°\-2R V'\x) (4.9) 

and we have used (3.21) to simplify some of these expressions. 

The coefficient c\, c 2 , C3 have a number of useful properties. First, from Eq.(3.24), it 
is easily seen that 

d = - 2a2.Ro 2 - 2b 2 4 (4.10) 

and thus c\ < 0. Second, using Eq.(3.23), 

c 2 = -2°\ (a 2 ai + b 2 R 2 ~) - tfb 2 °\ (4.11) 

Using (3.21), twice, it then follows that 

c 2 = -2a 2 a 2 x a 2 v - 2b 2 af > = c x (4.12) 

Third, c\ and C3 are related as follows. From Eq.(3.25), C3 may be written, 

( o 2 \ R 2 
c 3 = ARl a 2 + 4>b 2 = -2^Vi(4.13) 

V ° x J a * a P 
using (3.21) and (4.10). It follows that the linear terms in (4.6) may now be written, 

2R 2 

ciX + c 2 Y + c 3 Z = ci\X + Y - -^\Z J (4.14) 

a x a p 

Clearly (4.6) is zero at the stationary solution, but it cannot be negative for arbitrary 
X, Y and Z, because of the presence of the linear terms. However, X, Y and Z are not 
arbitrary but must respect the uncertainty principle (an expression of which is Eq.(3.19), 
for example). A convenient way to implement this restriction is to note that 

2R 2 

< (AA) 2 = a 2 (X + Y) J±Z (4.15) 

a x 
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with equality if and only if the state is a general Gaussian, such as the stationary solution. 
From (4.14), it is clear that 

Cl X + c 2 Y + c 3 Z = %(AA) 2 (4.16) 

Since c\ < 0, the linear terms are negative definite and zero only at the stationary solution. 
With some rearrangement of the quadratic terms, and using (3.21), 

M d -^ =^(AAf - - 2aV (X - Zf 



dt 




We therefore deduce that 



M d -^l < (4.18) 
dt ~ v ' 

with equality if and only if the solution is the stationary solution. This completes the proof 

of localization. 

As stated earlier, the stationary solutions to the Ito equation are valid for general 
potentials as long as the localization width is much less than the lengthscale on which 
the potential varies, i.e., as long as the approximation (2.20) holds. This approximation 
becomes exact for linear systems. 

We have essentially assumed the approximation (2.20) in proving the above localization 
theorem. This means that the proof is strictly valid only for systems with quadratic 
potentials. It cannot be valid for general potentials because even if there exist approximate 
stationary solutions for which the neglect of the higher derivative terms of the potential is 
valid, there will always be initial states for which (2.20) is not a valid approximation and 
localization is therefore not guaranteed for these states. For general potentials, therefore, 
the above proof implies localization only for a rather limited class of initial states, e.g., for 
states that are already close to the stationary states. 

Still, one intuitively expects that when approximate stationary solutions exist for gen- 
eral potentials, there will be situations in which most initial states will tend towards one 
of those solutions. Consider, for example, the case of a double well potential with minima 
a distance L apart, and suppose that the initial state has a width greater than L, where L 
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is chosen so that the approximation (2.20) is not valid. Then one can see from Eq.(2.11) 
that a very large initial width will be reduced very rapidly, in the mean, bringing it into 
the regime in which the approximation (2.20) is valid. Our localization theorem would 
then apply. We hope to investigate this point further in a future publication. 

Note that the stationary solutions and the localization theorem do not depend on the 
sign of V"(x), and therefore will be valid for the upside-down harmonic oscillator (which 
is sometimes used as a prototype for chaotic systems [41]). 



5. LOCALIZATION RATE 



It is also possible to estimate the rate of localization. Clearly, 

M dJAAf £ ci a 
at (jp 

and thus localization proceeds on a timescale of order r = a 2 /\c\\. Using (4.12), this 
becomes 

r = (2a 2 a 2 x + 2b 2 a 2 ^) (5.2) 

In the quantum Brownian motion model for the free particle with 6 = 0, Eqs.(2.5), (3.27) 
imply that 

r~(zLY (5-3) 



^kT J 

This, as noted previously, is the timescale on which thermal fluctuations become compa- 
rable to the quantum ones [42,43,44]. 

The above represents the minimum rate of localization. The actual rate can be much 
higher, e.g., if X is very large. Consider again the free particle with 6 = 0. Suppose, 
the initial state consists of a superposition of wavepackets a large distance £ apart. Then 

(Ax) 2 ~ e 2 , 

(AA) 2 pa Ah 2 \(3\ 2 (Ax) 2 ~ ^-f- (5.4) 
and the dominant contribution to the localization rate is the X 2 term, 

,rf(A^i) 2 „ 2 /„ 2 fi 2 \ „ 2 nW 
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It follows that in this case, 



For the quantum Brownian motion model, Eq.(2.5) then implies that 

%2 (5 7) 

firwykT 

Both of the timescales (5.3) and (5.7) are typically exceedingly small for macrosopic values 
of m, 7 and T. 

As we shall show in detail in the next section, once the solutions to the Ito equation 
have become localized, the corresponding density operator has the form (1.5). The local- 
ization timescale is therefore the timescale on which the density operator approaches the 
form (1.5). Since the process of decoherence of density operators is commonly associated 
with the approach to approximately diagonal form, it is natural to regard the localization 
timescale as essentially the same thing as the decoherence timescale. 

Note, however, that the so-called "decoherence timescale" is sometimes taken to be 
(5.7) [14,45,15]. What is clear from the above is that the rate of approach to diagonal form 
depends on initial state, and that (5.7) is appropriate only for initial states with very large 
(Ax) 2 . 

The connections between the timescales of decoherence and thermal fluctuations has 
certainly been noted before [43,42], but what is new here is the observation that both of 
these things are in turn related to the timescale of localization in quantum state diffusion. 



6. RECOVERY OF THE DENSITY OPERATOR 



We now show how a density operator satisfying the master equation may be recovered 
from the stationary solutions to the Ito equation. 

Each solution to the Ito equation is in general a functional of the noise term £(£) over 
the entire history of the solution's evolution. Eq.(1.2) indicates that the density operator 
is formally recovered from these solutions by averaging over all possible histories of 

the noise £(£), and we write 

p = M|^)(^| (6.1) 
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A completely explicit form of this expression may be found in Ref. [32] but it will not be 
needed here. 

When the solutions are the stationary solutions, (3.17), they depend on the noise 
£(£) only through their centres, (x), (p), which obey the Langevin equations (2.9), (2.10). 
We may therefore rewrite (6.1) as 

p = M J dpdq 8(p-p) 8(q-x) \ip pq ){ip pq \ (6.2) 

where we have again introduced the notation x = (x), p = (p), and \ip pq ) denotes the 
stationary solution (3.17) with centres p and q. The £(£) dependence is now contained 
entirely in p and x, and Eq.(6.2) may be trivially rewritten, 

p = J dpdq f{p,q,t) \ip pq ) (ip pq \ (6.3) 

where 

f(p,q,t) = M 8(p-p) 8(q-x) (6.4) 
The weight f(p, q, t) is non- negative and satisfies 

dpdq f{p,q,t) = 1 (6.5) 



It is in fact the solution to the Fokker-Planck equation corresponding to the Langevin 
equations. This Fokker-Planck equation is readily derived as follows. First note that 

f + df = M 8{p — p — dp) 8(q — x — dx) (6.6) 

Now expanding the delta functions to second order, we have 

f + df =M (5(p — p) 5(q — x) — dx 5(p — p) 5 f (q — x) — dp 5'(p — p) 5(q — x) 
H — dx 2 5(p — p) 8"(q — x) + dpdx 8'{p — p) 8' '(q — x) 

+±dp 2 S"(p-p)S(q-x)^ (6.7) 

We may now use the Langevin equations for x and p, and also pull the derivatives outside 
the mean, M, for example, 

M (dx 8(p - p) 8'(q -x))=M (j^KP - P) ^(Q - x)^) <& 

= ^ dt (6.8) 
m oq 
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We thus obtain the Fokker-Planck equation, 

df V df Tr „ ,df „ ,df 
at m oq op op 

+ k(P» L )| 2 ^2" + ^)| 2 ^i + 2Re (*(z, L)a(L,p)) ^j- (6.9) 

The coefficients of the second derivative terms are 

2 2 

|<r(p, L) I 2 = a 2 ^ + ft 2 ^ _ ft a & a 2 + ^ 

|(t(x, L)| 2 = 6 2 ^ + a 2 ^ - fta&a 2 + (6.11) 
2Re (a(x, L)a(L,p)) = 2a 2 J R a 2 + 26 2 a 2 - 2na6.Ro (6-12) 

We have 2%ab = 27, and for high temperature, the dominant term of the three second 
derivative terms is the first one, which has coefficient, 

\a(p,L)\ 2 « 2m-fkT (6.13) 

The resulting Fokker-Planck equation is well-known [46]. All solutions (for potentials for 
which e ~ v l kT is normalizable) tend towards the stationary solution, 

like e~^, where A is a normalization factor. For simplicity consider now the harmonic 
oscillator case V(q) = \muP'q 2 . Then the integrals over p and q may be done explicitly, 
with the result, 

/>(*, V) = exp f-!^(„ - yf - n ^yj-f \ ^- + y 2 )) (6.15) 



A v 37 2fcTA 
up to a normalization factor, where 

A = ^ + /3 + /3* (6.16) 

For large temperature, this is readily shown to be a thermal state [47] . Similar results are 
expected to hold for the case of more general potentials. 

To summarize, an initial density operator approaches the form (6.3) on the localization 
timescale, i.e., typically very quickly. On much longer timescales, it will then relax to an 
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equilibrium density operator, when one exists for the system (it does not for the free 
particle, for example). 

Note that although the above derivation of the asymptotic form (6.3) strictly concerned 
pure initial states, it is readily extended to mixed initial states by writing the initial state 
in a diagonal basis, 

P0 = Y^ °n \ n )( n \ ( 6 - 17 ) 
n 

and then applying the above to each initial state \n){n\. One thus finds that the density 
operator tends to the form (6.3), with f(p, q, t) of the form 

f(P, Q,t) = ^2 c n f n (p, q, t) (6.18) 
n 

where f n (p, q, t) is the solution to the Fokker-Planck equation corresponding to the initial 
state \n)(n\. 

As a final comment, note that any density operator may be written in the form (6.3), 
for some function f(p,q) - this is a property of the coherent states [36]. What is special 
about the particular function f(p,q,t) derived here is that it is non-negative, and that it 
obeys the Fokker-Planck equation (6.9). It may therefore reasonably be interpreted as a 
phase space probability distribution. (See Ref.[48] for related work on this point.) 



7. CONNECTION WITH THE 
DECOHERENT HISTORIES APPROACH 



As shown in Ref.[32], there is a close connection between the quantum state diffusion 
approach to open systems and the decoherent histories approach. In this section, we use 
the results of the previous sections to exemplify and amplify this connection. 

The primary mathematical aim of the decoherent histories approach is to assign prob- 
abilities to the possible histories of a closed system [4,29,30,31,49]. The approach is, 
however, applicable to open systems since they may be regarded as subsystems of a large 
closed system. A quantum-mechanical history is defined by an initial state pq at time 
t = to together with a string of projection operators P ai • • ■ Pa n acting at times t\ . . .t n , 
characterizing the possible alternatives of the system at those times. The projections are 
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exhaustive, ^2 a P a = 1 ; and exclusive, P a Pp = S aj 3P a - Due to interference, most sets of 
histories for a closed system cannot be assigned probabilities. The interference between 
pairs of histories in a set is measured by the so-called decoherence functional, 

D(a,a!) = Tr (P an (tn) • ■ ■ P ai (h) pP a [(h) • • ^<(*n)) (7.1) 

where P ak {^k) = e^^ Htk P a e^ Htk , is the Hamiltonian of the closed system and a denotes 
the string a\ • • ■ a n . When 

D(a,a')»° (7.2) 

for all pairs «/a', inteference may be neglected, and the set of histories is then said to 
be decoherent. One may then assign the probability p(a) = D(a,a) to the history, which 
may be shown to obey the sum rules of probability theory. 

For a given Hamiltonian and initial state, one's initial aim is to find those histories 
for which the decoherence condition is satisfied. In general, it is satisfied only by histo- 
ries which are coarse-grained, which loosely speaking, means that the projections at each 
moment of time give a less than complete description of the system. For open systems, a 
natural coarse-graining is to focus only on the properties of the distinguished system itself, 
whilst ignoring the environment. This involves using projections of the form, P a (g) I £ 
at each moment of time, where P a is a projection onto the distinguished subsystem and 
1^ denotes the identity on the environment. Assuming that the initial density operator 
factorizes, the trace over the environment may be carried out explicitly in the decoherence 
functional (7.1), and, in the regime in which a Markovian approximation holds, it then has 
the form 

D(a, cl) = Tt {PaK\U t^-i ' ' ' < 2 [^{j [Po]^J ' ' " P a' n JPa n ) (7-3) 

where the trace is now over the distinguished subsystem only. The quantity is 
the reduced density operator propagator associated with the master equation (1.1), pt = 

The results of the previous sections have provided us with some information about the 
density operator propagator, and we can use this information to establish some properties 
of the decoherence functional (7.3). 

We have seen that any density operator will tend, on a typically very short timescale, 
to the form (6.3), in which it is approximately diagonal in a set of phase space localized 
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states. Once in that form, under further evolution its form will be preserved and the 
only change will be that the function f(p 7 q, t) will evolve according to the Fokker-Planck 
equation (6.9). 

Take the projection operators in the decoherence functional to be phase space projec- 
tors, of the form 

P a = dqdq\i) pq )(i) pq \ (7.4) 

where \ip pq ) are the generalized coherent states (3.17), and are eigenstates of the operator 
(4.1). These quantities are not exact projection operators, but will be approximate pro- 
jectors if the phase space region T a is sufficiently large, and if its boundary is sufficiently 
smooth [30]. They have the property that P a \^pq) ~ \ippq) if P,Q he in the phase space 
cell r a , and P a \i^ P q) ~ otherwise. Again this approximation should be valid if r Q is 
sufficiently large compared to the phase space area occupied by the generalized coherent 
states (which is of order h) . 

Consider the time evolution from to to t\ in the decoherence functional. Clearly if this 
time interval is greater than the localization time it follows from the results of Section 6 
that the density operator will evolve into the form 

K tl[po\ = J d P d <l fiPi^hMpqtiippql (7.5) 

Because it is approximately diagonal in the coherent states, it is easy to see that 

P ai ^[p]P a; «0 (7.6) 

if ol\ 7^ a^. Keeping only the diagonal terms, ot\ = aj, and evolving to time t^-, the 
(unnormalized) density operator P ai K^[p]P ai should again evolve into approximately 
diagonal form, and again we get 

P a > 2 « (7.7) 

if «2 7^ a^. Continuing in this way for the entire history, it is easy to see that we will have 
approximate decoherence if the projections at each moment of time are taken to be phase 
space projectors. We have not estimated the degree of approximate decoherence (and this 
tends to be rather involved in general), but we expect it to be good if the size of the 
phase space cells is much larger than %, and if the time between projections is longer than 
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P K t2 



P ai K\l\p\P ax 



the localization time. We therefore find that localization in quantum state diffusion and 
decoherence of histories in the decoherent histories approach occur in the same variables. 

This conclusion is in agreement with the general connection between quantum state 
diffusion and decoherent histories outlined in Ref.[32], but it also extends it somewhat. 
There, it was argued that localization and decoherence tend to occur in the Lindblad op- 
erators. Here, the Lindblad operator is essentially position, but we have actually obtained 
the stronger conclusion that localization/decoherence occurs in the operator (4.1), and 
hence, approximately, in both position and momentum. (Note that the Lindblad operator 
has a small momentum part added, but this is not the primary source of momentum lo- 
calization. Rather, it is the interplay between the position part of the Lindblad operator 
and the Hamiltonian, as discussed earlier). 

Given approximate decoherence, we now consider the probabilities for histories, given 
by the diagonal elements of the decoherence functional. From Eq.(7.5), and from the 
properties of the phase space projections, it follows that 

p aiKttM p ai ~ / dpidqi f(pi,qi,h) \^ Pm ) (^ Piqi \ (7.9) 

Jr ai 

Now consider the evolution from t\ to t<i- We have, from Section 6, 

K tl [l^PiqJ&Piqil] = J d P2dq2 f(P2, 92, ^bl, 91, *l) iVfe^KV^I ( 7 ' 10 ) 

where f(p2, 92, *2 |pi > 91, ^l) is the solution to the Fokker-Planck equation satisfying the 
initial condition, 

/(P2»S2»*lbl»?l»*l) = KP2-Pl) $(42 - 9l) (7-11) 

/(P2, 92, ^2 |Pl 5 91) ^l) is therefore the Fokker-Planck propagator, i.e., the probability of 
finding the particle at P2, 92 at time £2, given that it was at pi, q\ at time t\. Assembling 
(7.9) and (7.10), it follows that 

Pa 2 K tl P^K^ipolP^ Pa 2 ~ dp 2 dq 2 / dpidqi f(p2, 92, t 2 \pi, 91,^1 ) 

X /(Pl,9l,*l) \^Piqi)(^Piqi\ (7-12) 

Continuing in this way for the entire history, one finds that 

p(ai, • • ■ ctn) = dp n dq n ■■■ / dpidqi f{p n , 9n, t n \p n -l, 9n-l, *n-l) 

"'Tan J^ai 

X ••• /(P2,92,*2bl,9l,*l) /(Pl,9l,*l) (7-13) 
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This is the desired result. Eq.(7.13) is the probability that a particle evolving according 
to the stochastic process described by the Fokker-Planck equation (6.9) will be in the 
sequence of phase space cells r ai • • • T an at times t\- ■ -t n . 

This result is in agreement with the probabilities one would assign to histories in the 
quantum state diffusion approach. For there, once the solutions to the Ito equation have 
become localized, the description of the motion on scales greater than the localization 
width is classical Brownian motion according to the Langevin equations (2.9), (2.10). This 
is equivalent to the description in terms of the Fokker-Planck equation (6.9). We have 
therefore exemplified the second part of the connection between quantum state diffusion 
and decoherent histories put forward in Ref. [32] - that the probabilities assigned to histories 
in each approach are the same. 

A further claim in Ref. [32] is that the degree of localization is related to the degree of 
decoherence. Although they are clearly related, it is difficult to check this here because, 
as stated above, explicit computation of the degree of approximate decoherence is quite 
difficult. This point will be pursued in more detail elsewhere. 

Finally, a property of the Fokker-Planck propagator associated with Eq.(6.9) is that it 
is peaked about classical evolution (with dissipation). It follows that the probability for 
histories (7.13) will be most strongly peaked when the phase space cells lie along a classical 
path. 

8. SUMMARY AND DISCUSSION 

Our main results follows. 

We have found stationary solutions to the Langevin-Ito equation (2.4) which are exact 
for linear systems, and approximate for non-linear systems as long as the localization width 
is much less than the scale on which the potential varies. The solutions consist of localized 
wave packets concentrated about a point in phase space undergoing classical Brownian 
motion. 

For linear systems, every initial state tends towards one of the stationary solutions. 
For non-linear systems, some form of localization is plausible, and will certainly be true 
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in the neighbourhood of the stationary solutions, but our investigations on this point are 
inconclusive. 

Localization proceeds on a timescale which is typically very short. It is related to the 
timescale on which thermal and quantum fluctuations become comparable, and also to the 
decoherence timescale. 

The density operator corresponding to the stationary solutions may be reconstructed 
and has the form (1.5). It is therefore diagonal on a set of phase space localized states. 
For linear systems (and plausibly for many non-linear systems also) any initial density 
operator approaches this form on the localization time scale. On longer timescales, when 
dissipation is present, the density operator approaches a thermal state (when it exists) in 
the long-time limit, as expected on general grounds. These results fulfil the aims set out 
in the Introduction, concerning the density matrix approach to decoherence. 

Our work also has some implications for the question of approximate versus exact 
density matrix diagonalization. As discussed in the Introduction, it is often held important 
in the context of decoherence studies to find the basis in which the density matrix is 
diagonal. This can of course always be done, since the density operator is a hermitian 
operator, but the basis in which p is exactly diagonal is generally non-trivial, i.e., it does 
not usually consist of the eigenstates of a simple operator. Furthermore, the basis consists 
of eigenstates of a different operator at each moment of time. 

Here, we have shown that the quantum state diffusion approach naturally leads to a 
basis in which the density matrix is approximately diagonal. The basis states are the eigen- 
states of a simple operator, the same operator at each moment of time. There therefore 
appears to be much to be gained by relaxing the condition of exact diagonality. Corre- 
sponding to these exactly and approximately diagonalizing bases, there will be exactly and 
approximately decoherent set of histories in the decoherent histories approach. In Section 
7, we exhibited the approximately decoherent set. 

The bases of approximate and exact diagonality do not appear to be "close" in any 
sense. For example, for a Gaussian density operator (in the position representation), the 
exactly diagonal basis consists of Her mite polynomials multiplied by Gaussians (similar 
to energy eigenstates of the harmonic oscillator) [5], whereas the approximately diagonal 
one consists of phase space localized states. (See also Ref.[50] for examples of different 
bases in which the density matrix is diagonal). This suggests that the corresponding 
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exactly decoherent set of histories is not necessarily "close" to the approximately diagonal 
one, somewhat contrary to the expectation sometimes expressed [51] (although it is not 
clear whether there are other exactly decoherent sets of histories that are close to the 
approximate one). 

The basis of states picked out by the QSD approach appears to be "natural", in the 
sense that they correspond to the trajectories that would actually be observed in an individ- 
ual experiment, whereas the exactly diagonal basis does not, in general. Correspondingly, 
the approximately decoherent set of histories may seem to be more "natural" than the 
exactly decoherent set. The question of whether one is any sense preferred over the other 
is, however, a subtle one. It depends on the sort of predictions one wishes to make, and on 
the extent to which the simplified situation consisting of a distinguished system coupled 
to an environment is really part of a much larger universe in which there may be adaptive 
systems that can measure different properties of the distinguished subsystem [4] . 

The sum up, the model described in this paper illustrates the connection between the 
intuitive pictures and physical predictions provided by the quantum state diffusion ap- 
proach, density matrix approaches, and the decoherent histories approach. In our model, 
localization in quantum state diffusion, diagonalization in the density matrix approach, 
and decoherence of histories in the decoherent histories approach all occur under the same 
conditions and are essentially the same thing, for each is concerned with the conditions 
under which "definite properties" may be assigned to the system. Furthermore, the prob- 
abilities assigned to histories in the quantum state diffusion approach and the decoherent 
histories approach approximately coincide. 
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